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ABSTRACT 

We present our spectroscopic data around H/? of the symbiotic star V1016 Cyg 
obtained with the Bohyunsan EcheUe Spectrograph, in order to secure the broad 
emission feature at around 4850 A, which is formed through Raman scattering of 
He II A 972. The total cross section around Lyj is approximately given by a{\) ~ 
1.7 X lO^^^lX^yy/ (A — Aiy-y)]^ cm^, with Xlvj being the line centre wavelength of Lyj. 
We find a centre shift redward by an amount AA = +0.64 A in the Raman scattered 
He II A 4850. This redward centre shift is exactly analogous to the effect for the Raman 
scattered He II A 6545 blueward of Ha discussed in our previous study. We compute 
the branching ratios of Raman scattering into the level 2s and levels 3s, and 3d, which 
are subsequently incorporated in our Monte Carlo code. Using this code, we present 
the centre shift of the 4850 feature as a function of the neutral hydrogen column 
density in the scattering region. Assuming that He II A 972 emission is characterized 
by a single Gaussian profile, the redward peak shift observed in the Raman scattered 
He II 4850 feature corresponds to the neutral column density Nhi = 1.2 x 10^^ cm~^. 
Assuming that the covering factor ~ 0.1 of the scattering region with respect to the 
He II emission region and adopting a simple spherical stellar wind model, we may 
place an upper bound M ^ 3.6 x 10~^ Mq yr~^ for the mass loss rate of the giant 
component of V1016 Cyg. Our estimate can be severely affected by the kinematics of 
the scattering and He II emission regions and the exact atomic physics, about which 
brief discussions are presented. 

Key words: scattering — radiative transfer — binaries: symbiotic — mass-loss — 
individual V1016 Cyg 



1 INTRODUCTION 

Symbiotic stars, believed to be a binary system of a giant 
and a hot white dwarf, are important objects for studying 
mass loss processes occurring in the late stage of stellar evo- 
lution (e.g. Kenyon 1986). A thick neutral hydrogen region 
formed around the giant component as a result of the mass 
loss is illuminated by strong UV radiation originating from 
the hot white dwarf component. Schmid (1989) identified the 
broad emission features around 6825 A and 7088 A found in 
more than half of the symbiotic stars by suggesting that they 
are formed via Raman scattering of the resonance doublet 
O VI 1032, 1038 by atomic hydrogen. Simultaneous obser- 
vations in the far UV regions and many spectropolarimetric 
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observations strongly support this identification (e.g. Bir- 
riel, Espey & Schulte-Ladbeck 1998, 2000, Schmid & Schild 
1994, Harries & Howarth 1996). 

Raman scattering also appears to operate in the forma- 
tion of Ha wings that have been observed in most symbiotic 
stars and in a number of young planetary nebulae (Nuss- 
baumer, Schmid & Vogel 1989, Lee 2000, Lee & Hyung 2000, 
Arrieta & Torres- Peimbert 2003) . The Ha profiles in the far 
wing regions are well fitted by fx oc AA~^ = (A — XHa)~^, 
which delineates approximately the Raman scattering cross 
section, where Xho is the line centre wavelength of Ha. He II 
n ^ 2 emission lines arising from states with even principal 
quantum numbers have wavelengths only a little bit shorter 
than those of the H I Lyman series lines. These He II emis- 
sion lines can be Raman scattered to form broad features 
blueward of the Balmer series lines. 

We expect that the strongest Raman scattered He II 
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feature is formed at around 6545 A blueward of Ha, which is 
formed through Raman scattering of He H A 1025. This fea- 
ture was found in the spectra of the symbiotic stars RR Tel, 
He 2-106 and V1016 Cyg (Lcc, Kang & Byun 2001, Lee et 
al. 2003). Another broad feature is formed at around 4850 
A, which is the Raman scattered feature of He H A 972. This 
feature was reported in spectra of the symbiotic stars RR 
Tel (Van Groningen 1993) and V1016 Cyg (Birriel 2003). It 
is quite notable that the young planetary nebula NGC 7027 
also exhibits this 4850 broad feature (Poquignot ot al. 1997). 

Raman scattered He II lines can be used as a powerful 
diagnostic of the mass loss processes in symbiotic stars, be- 
cause the mass loss rate can bo calculated by inferring the 
optical depth structures and the covering factor of the scat- 
tering region with respect to the UV emission region around 
the white dwaxf. However, the flux of a Raman scattered 
He II feature is mainly determined by the product of the 
optical depth and the covering factor, which is a big barrier 
for the exact estimate of the mass loss rate. 

In the previous theoretical study by Jung & Lee (2004, 
hereafter JL04), they showed that the Raman scattered 
He II A 6545 should exhibit a redward centre shift by a 
significant amount depending on the H I column density 
of the scattering region around the giant in the absence of 
the relative motion between the scattering region and the 
emission region, which will lift the degeneracy of the opti- 
cal depth and the covering factor. However, Raman scat- 
tered He II 6545 is usually blended with the forbidden line 
[N II] A 6548, making it very hard to determine the ex- 
act peak location. This difficulty may be avoided when we 
use the Raman scattered He II 4850, which should be much 
weaker than the 6545 feature but be much easier to locate 
the peak wavelength. 

UV radiation around Ly7 may have four scattering 
channels depending on the final state of the scattering 
atom, which can be one of Is, 2s, 3s and 3d. When scat- 
tering occurs into either 3s or 'id level, energy conser- 
vation leads to the wavelength of the scattered radiation 
Ass.sd = (A^g ji - ^lIp)'^ = 1-85 M, where \hr ii 972 
and Xhyii are line centre wavelengths for He II A 972 and 
Ly/3. Therefore, in order to simulate the Raman scattered 
He II lines accurately, it is necessary to compute the branch- 
ing ratios into 3s and id states. This calculation will also be 
useful for the future study of Raman scattered lines formed 
blueward of Paschen series lines. 

In this paper, we present an optical spectrum around 
H/3 of the symbiotic star V1016 Cygni obtained with the 
1.8 m Bohyunsan telescope and our Monte Carlo calcula- 
tions of Raman scattering around Ly7, in order to deter- 
mine the neutral hydrogen column density of the scattering 
region around the giant component. 



2 OBSERVATIONAL DATA 

We observed V1016 Cyg on the night of 2003 December 16 
with the Bohyunsan Optical Echelle Spectrograph equipped 
with a 2kx4k E2V CCD installed on the 1.8 m telescope at 
Mt. Bohyun. We used the 200 ji optical fiber that yields the 
spectroscopic resolution R = 44, 000 and the exposure time 
was 3000 s. The data have been reduced following standard 
procedures using the IRAF packages. In Fig. 1, we show the 
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Figure 1. The optical spectrum around H/3 of V1016 Cyg ob- 
tained with the Bohyunsan Optical Echelle Spectrograph (BOES) 
installed on the 1.8 m telescope at Mt. Bohyun. The spectroscopic 
resolution is 44,000 and the exposure time was 3000 s. The thin 
soUd lines represent our single Gaussian fits to the emission Unes 
H/3, He II A 4861 and the Raman scattered He II A 4850. From our 
Gaussian fitting, the Raman scattered He II A 4850 is shifted red- 
ward by an amount AA = -1-0.64 A with respect to the emission 
fine center He II A 4861. 

part of our spectrum that contains H/3. The solid line shows 
the observed spectrum, and in the inset we enhanced the 
part exhibiting the Raman scattered He II A 4850 feature. 
In this figure, the vertical a^xis is shown in an arbitrary unit, 
which is sufficient for the purpose of this work. 

We measure the relative shift from the emission lines 
He II A 4861 and H/3, where these emission lines are assumed 
to be isotropic in their emission region. This assumption may 
not be true in the accuracy required in this work. In the case 
of O VI A 1032, 1038 doublet, it is known that their profiles 
and those of their Raman scattered 6825, 7088 features differ 
quite significantly (Schmid et al. 1999). The O VI doublet 
lines arc strong resonance lines that are subject to various ra- 
diative transfer efi^ects, whereas optical He II recombination 
lines are expected to be relatively free from such effects. This 
may be checked by monitoring the variations as a function 
of the orbital phases. In this work, we take the He II A 4861 
and H/3 as the local reference of velocity and measure the 
line centre shift of the Raman scattered He II A 4850. 

In Fig. 1, we applied single Gaussian fits to both 
He II A 4861 and H/3 to locate their peak positions. The 
dotted lines show our fit, which we can barely notice in 
the figure. This implies that satisfactory fits are obtained 
for both emission lines. For He II A 4861, the fitting func- 
tion is /we // = /le-K^-^i^/'^^il' and for H/3 it is fnn = 
f^(.-[{>^-^2)/£^>~2? ^ where /i = 165000, /g = 8200, Ai = 
4860.54 A, A2 = 4858.57 A and AAi = AA2 = 0.75 A. 
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merit for bound-bound state transition is given by 



^ = in'/6 (pole) 
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Figure 2. The contour C for the integral in Eq. 9. There is a 
branching cut on the real axis between ^ = 1/2 and £, = —1/2. 
There is a pole at ^ = in' 1%, at which residues axe computed. 



The 4850 feature is also fitted with a single Gaussian 



he 



-[(A-A3)/AA3]^ 



(1) 



where the least squaxe method yields As = 4850.63 A and 
AA3 = 3.76 A. We note that the width broadening is in 
excellent agreement with the theoretical expectation 



AA3/AA1 ~ 4861/972 = 5.0. 



(2) 



Using the atomic data compiled by van Hoof, the centre 
wavelength Ao in vacuum of the Raman scattered He H A 972 
is 



Ao = (A^^ 



- A^^)"^ = 4852.098 A, 



(3) 



where the vacuum wavelength of Lya is Aa = 1215.67 A and 
that for He II A 972 is Xne 11 972 = 972.112 A. Considering 
the refractive index of air Uair = 1.000279348, the line cen- 
tre of the 4850 feature appears at A3 = 4850.743 A in air. 
Noting that the line center of He II A 4861 in our spectrum 
appears at \h&ii = 4859.320 A, our spectrum of V1016 Cyg 
shows the redward centre shift AA = -1-0.64 A in the 4850 
feature with respect to their emission region. 



3 ATOMIC PHYSICS 

We compute the matrix elements necessary for computation 
of the contribution to the scattering cross section of the 
levels 3s and 3d. The contributions of the levels Ls and 2s 
have been presented in the literature (e.g. Saslow & Mills 
1967, Sadegpour & Dalgarno 1992, Lee & Lee 1997). 

According to Karzas & Latter (1961), the matrix ele- 
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n — n 
n + n' 



-Ann' /{n — n')^ and F(q,/3, 7, a:) is the hy- 



pergeometric function 

Here, Rni (r) is the radial wavefunction of a hydrogen atom, 
which can be expressed using a confluent geometric function 
Fi by 



Rniir) 



(2/-I-1)! V {n-l-iy.2n \ n J 
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where Fi{a,b;x) = 1 + abx + [a{a+ 1)&(6-1- l)/2!]x^ -I . 

By setting n' = 3, Z = 1, a straightforward calculation 
leads to 
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By interchanging n and n' in Eq. (1) and letting n' 
3,1 = 2, we obtain 
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According to Bethe & Salpeter (1957) the wavefunction 
for the continuum states n'p is given by the contour integral 
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4 Jung & Lee 



L 
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where the contour C is depicted in Fig. 2. 

The matrix element for a dipole operator between 3s 
and n'p states is 
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Using the residue theorem for the pole at ^ = in /6, we 
compute the contour integral 
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(n'2 + 9)4 
Substituting this result we have 

^ 3^/^12^n'^/^(l + n'^)^/^(7n'^+27) _2„'tan-(3/n') 

^3. (l_e-2Tn')(„/2+9)4 e 

A similar procedure gives 
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In Fig. 3 wc show the scattering cross sections and 
branching ratios as a function of the wavelength of UV ra- 
diation around Ly7. In the upper panel we show the total 
scattering cross section ax by a thick solid line, the cross 
section (T2s for Raman scattering into the 2s level by a dot- 
dashed line, and the cross section a^sSd for Raman scattering 
into the 3s — 3d levels by a dashed line. 

Also in the upper panel, the dotted line shows a 
Lorentzian fit to the total scattering cross section, which 
is given by 



a (A) = (Jo 



A- A 



(15) 



Here, Alj/7 is the line centre wavelength for the Ly7 transi- 
tion. With this fit, we set ao = 0.17x lO'^'^cm^. As is for the 
case of Lya shown by Lee (2003), the total cross section is 
also asymmetrical with respect to the Ly7 line center, which 
is noticeable in this fit. At the line centre A = 972.112 A of 
He II A 972, the total scattering cross section is 



(Ttot.He II 972 = 0-91 X 10 Cm^. 



(16) 



This may be compared with the total scattering cross section 
for He II A 1025, which is <Jtot,He ii 1025 = 6.2 x 10~^^ cm"^. 
From this computation, we may regard the Raman scattered 
He II A 6545 as a probe for neutral regions with Nhi ~ 
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Figure 3. The scattering cross sections and branching ratios 
of UV radiation around Ly7 by atomic hydrogen. In the upper 
panel, the solid line shows the total scattering cross section, the 
dot-dashed line represents the scattering cross section into the 
level 2s (for He II A 4850) and the dashed line represents the 
sum of the scattering cross sections into the levels 3s and 'id. 
Also in the upper panel, the dotted line shows a Lorentzian fit 
to the total scattering cross section by the formula ctq [^Lyi/ (-^ ~ 
^Ly-i\~^-i where \Lyy is the line centre wavelength of Ly7 and 
(To = 1.7 X 10~2* cm2. In the lower panel, we show the branching 
ratio ri of Raman scattering into the level 2s by a solid line, and 
the branching ratio r2 into the levels 3s and 3d by a dashed line. 



10^" cm~2 and the Raman scattered He II 4850 is formed 
mainly in regions with Nhi ^ 10 cm~ . 

In the lower panel of Fig. 3, we show the branching ra- 
tios ri = a2s/(7r,f2 = aSsSd/o'T for Raman scattering into 
the 2s level and the 3s— 3d levels by a solid line and a dashed 
line, respectively. It is quite notable that the branching ratio 
into the level 2s is about 3 times larger than the sum of the 
branching ratios into 3s and 3d near He II A 972. Raman 
scattering into levels 3s and 3d will leave a broad feature at 
A = 1.85/i blueward of Paa. According to this calculation, it 
is expected that the A = 1.85/i feature will be constituted by 
a smaller number of photons by a factor of 1/3 than those 
in the Raman scattered He II A 4850. At around A = 985 A 
the total cross section is dominated by the scattering chan- 
nel into the levels n = 3, for which, however, the scattering 
cross section is very small. Near the line center of Ly/3 the 
branching ratio for the levels n = 3 becomes negligibly small 
due to the decrease of the phase space volume and becomes 
zero redward of Ly/3, which is forbidden by energy conser- 
vation. 
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4 CALCULATIONS 

4.1 Monte Carlo Procedures 

The same Monte Carlo code used in the previous study by 
JL04 is slightly modified with an addition of subroutines for 
those scattering channels into 3s and 3d levels. The code in- 
corporates the channel for Rayleigh scattering as well as the 
two channels for Raman scattering. Therefore, we followed 
faithfully paths of He II photons suffering several Rayleigh 
scatterings followed by a Raman scattering in the H I region. 
However, no dust effect is considered in the current work. 

As in JL04, considering the fact that He II 6560, 
He II 4686 emission lines are well fitted by a single Gaussian 

fx = /oexp[-(A - Ao)VAA2], (17) 

with the velocity width AA/Ao = (23 km s~^/c). We use 
the same single Gaussian for the incident radiation. Since 
we are only concerned with the exact location of the central 
peak, we do not require the exact equivalent width of the 
He II A 972 line and set it to be O.lA. 

When the scattering region is optically thin, the 
strength of the Raman scattered feature is mainly deter- 
mined by the product of the incident radiation flux and the 
scattering cross section. The He II A 972 incident line radi- 
ation has a profile that is symmetric with respect to its line 
center, whereas the scattering cross section is highly inclined 
toward the line center of Ly7, which is located redward of 
He II 972. This induces a redward shift of the peak in the 
Raman scattered He II A 4850, as is explained for the Raman 
scattered He II A 6545 in JL04. We illustrate this situation 
in Fig. 4, whore the dotted line shows the Raman scattering 
cross section into the level 2s relevant to the formation of the 
4850 feature. In this figure, the dashed line shows the sum 
of the Raman scattering cross sections into the levels 3s and 
3d, which is slightly larger than the quantity represented by 
the dotted line in this region. 

Due to the multiple scattering effect that is important 
for a very optically thick scattering region, the peak shift 
will decrease as Nhi increases, and therefore a Monte Carlo 
technique is useful to obtain a reliable quantitative relation. 

4.2 Results 

4-2.1 Balmer and Paschen Wings 

In Fig. 5, we show typical Raman scattered He II features 
formed in the blue region of Ha, H/3 and Paa, where we set 
Nhi = 3 X 10^° cm"^. The formation of the 6545 feature 
was studied intensively by JL04, and here it was reproduced 
for comparisons with those obtained from Raman scattering 
of UV radiation around Ly7. The vertical axis shows the 
number of photons obtained from the Monte Caxlo simula- 
tions for 5 X 10'' continuum photons per wavelength interval 
of AA = 0.00125 A in the UV rest frame around Ly7. 

We normalize the horizontal scale by the velocity 
c(AA_r/A_r) from the lino centres of Raman scattering, where 
Ah = 6563,4861, and 18751 A respectively. With this nor- 
malization, the flux density at Xr in the Raman scattered 
features is diluted by the factor XR/Xhyy- 

Far from the line centres of Balmer lines and Paa, the 
Rayleigh and Raman scattering optical depth is quite small, 
and the Raman wings exhibit profiles that are proportional 
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Figure 4. The Raman scattering cross sections overplotted to 
the He II 972 emission line given by a single Gaussian with AA = 
0.076 A. The dotted line shows the cross section of the Raman 
scattering into the 2s level, and the dashed line shows the sum 

of the cross sections for the Raman scattering into the ,36' and 
Sci levels. Note that the cross section increases steeply toward the 
Ly7 line center, which will affect the center location of the Raman 
scattered He II 4850 feature. 

to the total scattering cross section multiplied by the corre- 
sponding branching ratio. It is notable that the sum of the 
branching ratios for levels 3s and 3d is smaller by a factor 
1/3 than the branching ratio for level 2s, and therefore our 
simulations show that the Raman scattered feature blue- 
ward of Paa is weaker by the similar factor than the feature 
blueward of H/3. 

Near the line centres of Balmer lines and Paa, the total 
scattering optical depth is quite large and the Raman con- 
version efficiency no longer increases toward the line center. 
Instead, almost all incident Ly7 photons may be singly scat- 
tered at the surface of the scattering region by a Rayleigh 
or Raman process, or may suffer a few Rayleigh scatterings 
followed by Raman scattering or final Rayleigh scattering 
before escape, which leads to nearly constant Raman con- 
version efficiencies into around H/3 and Paa. More quan- 
titatively, according to our Monte Carlo results, about 36 
percent of incident photons very near Ly7 line center are 
converted and redistributed around H/3 and 13 percent of 
them are redistributed around Paa. 

It is also quite notable that the Raman feature blue- 
ward of Paa appears 4 times further away from the centre 
of Paa than the 4850 feature is away from the H/J centre. 
This is attributed to the difference in the ratio Ao/Aj of the 
wavelengths of the Raman scattered radiation and incident 
UV radiation. 

Prom this figure, we infer that the wings around H/? and 
Paa are very weak compared with those around Ha, unless 
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Figure 5. Typical Raman scattered He II features formed blue- 
ward of Ha, up and Paa. The H I column density of the scat- 
tering region is assumed to be Nhi = 3 x 10^° cm~^. We nor- 
malize the horizontal scale by the velocity c{AX[{/\ii) from the 
line centres of Raman scattering, where = 6563,4861, and 
18751 A respectively. Note that with this normalization the flux 
density is diluted by the factor Xji/XLyy - The vertical axis shows 
the number of photons obtained from the Monte Carlo simula- 
tions for 5 X 10* continuum photons per wavelength interval of 
AA = 0.00125 A in the UV rest frame around Lyj. 



the continuum around Ly7 is unrealistically stronger than 
that around Ly/3. The relative strengths of Ha and H/3 wings 
may be useful in distinguishing those wings that are formed 
from other mechanisms including Thomson scattering. 



4-2.2 H/3 wings and the Raman Scattered He II 4850 

In Fig. 6 we show H/3 wings and the Raman scattered 
He II A 4850 feature for 4 values of Nhi- The extent of 
H/3 wings increases as Nhi increases and the wing profile 
is approximately by (A — \h/3)~^- A saturation behaviour is 
seen near the line center of H/3 at which the Raman scatter- 
ing optical depth becomes larger than unity. The strength 
of the Raman scattered He II X 4850 also increases as Nhi- 

In Fig. 7, in order to illustrate the centre shift phe- 
nomenon clearly, we overplot profiles of the Raman scat- 
tered He II A 4850 obtained for various Nhi as we in- 
crease Nhi starting from Nhi = 2 x 10^° cm~'^ by a step 
of ANhi — 4 X 10^" cm~^. Here, we use the same incident 
radiation as in Fig. 6. Because the strength of the 4850 fea- 
ture increases as Nhi, stronger profiles correspond to higher 
Nhi- It is immediately seen in this figure that the amount 
of centre shift decreases as Nhi- 

In Fig. 8, we show the centre shift of He II A 4850 fea- 
ture as a function of Nhi, treating the scattering region as 
a plane-parallel slab characterized by a single column den- 



Figure 6. The Raman scattered He II A 4850 feature formed in 
neutral scattering regions with various H I column densities. The 
strength and the centre shift of the He II A 4850 feature exhibit 
similar behaviours to those of the He II A 6545 feature discussed 
by JL04. The incident radiation is prepared in the same way as 
in Fig. 5 and the vertical axis has the same meaning as in Fig. 5. 



sity. As is emphasized earlier, we assume that the scattering 
region is static with respect to the emission region and mea- 
sure the relative shift from the He II A 4861. There exist 
fluctuations in the result shown in Fig. 8, which is indica- 
tive of the statistical nature of the Monte Carlo calculation. 
We provide a simple fit to our data, which is given by 



(AA/1 A) =0.73 



0.015a: + 0.83(0.8 (18) 
cm~'^ is the H I column density in 



-,20 



where x = Nhi/W' 
units of 10^° cm"^. 

The behaviour of the centre shift for He II A 4850 as a 
function of Nhi is almost similar to that for He II A 6545. 

For Nhi ^ 10^'' cm^^, the amount of centre shift decreases 
as Nhi increases. We omit the description of our results, 
which can be found for Raman scattered He HA 6545 in 
JL04. According to our Monte Carlo calculation, the ob- 
served centre siiift by an amount 0.64 A is obtained when 
Nhi — 1-2 X 10^^ cm~^. It is particularly notable that the 
slope of the centre shift as a function of Nhi is small for 
high column densities, so that the accurate determination 
of Nhi may be significantly affected by the relative motion 
between the He II emission region and the neutral scattering 
region. 



4.3 Estimate of the Mass Loss Rate 

Assuming a spherical stellar wind from the giant component 
in V1016 Cyg, we may set the hydrogen number density 
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Figure 7. The profiles of He II A 4850 feature obtained for 
various column densities overplotted. The incident He II A 972 
is assumed to be given by a single Gaussian with a width AA = 
0.076 A and the neutral scattering region is a plane-parallel slab 
characterized by a single column density Nhi- Nhi is increased 
from Nhi = 2 x 10^° cm'^ in steps of AAT^j = 4 x 10^° cm-^. 



Nj,, (10=0 cin-=) 

Figure 8. The centre shift of He II A 4850 feature as a function of 
Nfji- The incident He II A 972 is assumed to be given by a single 
Gaussian with a width AA = 0.076 A and the neutral scattering 
region is a plane-parallel slab characterized by a single column 
density Nfji- We give a simple fit to our data by (AA/1 A) = 
0.73 - 0.015a; + 0.83(0.8 + x)-°-^ , where x = Nhi/W'^° cm'^ is 
the H I column density in units of 10^" cm~^. 
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Figure 9. A schematic diagram for the scattering geometry in 
V1016 Cyg. Depending on the density structure around the giant 
and the strength of the UV radiation from the white dwarf, the 
Stromgren boundary is formed, which is defined as the loci of 
the ionization front dividing the H I region from H II region. 
The point H is the foot of the perpendicular to a siglitlinc from 
the white dwarf with an impact parameter p making an angle So 
with the axis connecting the giant and dwarf. The parameter s 
measures the distance from H to & given point along the sightline 



with a distance r from the giant so that 



■ p . By ds we 



denote the smallest angle of the line of sight from the white dwarf 
making with the axis connecting the giant and the white dwarf, 
for which no H I region can be reached. 



n(r) = 



M 



(19) 



where M is the mass loss rate of the giant, nip is the proton 
mass, n is the mean molecular weight, and Voo is the terminal 
velocity of the slow stellar wind. 

Depending on the density structure and the ionizing lu- 
minosity, the ionization structure will be determined (e.g. 
Taylor & Seaquist 1984). A schematic diagram for this situ- 
ation is shown in Fig. 9, where d is the binary separation. In 
this figure, we denote the ionization front by the Stromgren 
boundary which separates the H I region around the giant 
from the H II region formed around the white dwarf. By 6s 
we denote the smallest angle of the line of sight from the 
white dwarf mating with the axis connecting the giant and 
the white dwarf, for which the H I region cannot be reached. 

We consider a sightline from the white dwarf to the H I 
region. The point H is the foot of the perpendicular to this 
sightline from the white dwarf with an impact parameter p 
making an angle 6o with the axis connecting the giant and 
dwarf. The parameter s measures the distance from H to a. 
given point along the sightline with a distance r from the 
giant so that r'^ = s'^ + p^. The H I column density for this 
sightline with the impact parameter p is given by 



N(p) 



f 



n{r)ds 



M 



A'KiimpVooP 



TT _l SO 

- -I- tan — 
2 p 



(20) 



8 Jung & Lee 



where so is the distance from H to the Stromgren boundary. 

Brocksopp et al. (2002) analysed optical and radio im- 
ages of V1016 Cyg obtained with the Hubble Space Tele- 
scope and the Very Large Array (VLA). Adopting a distance 
of 2 kpc they suggest that the projected binary separation 
of V1016 Cyg is d = 84 ± 2 AU. Previous works by other 
researchers provide smaller binary separations and shorter 
binary periods. For example, Schmid & Schild (2002) per- 
formed spectropolarimetry to investigate the variation of the 
polarization direction which is expected to vary due to the 
orbital motion. Bcised on their polarimetric method, they 
proposed that the binary period may exceed 100 years. Lee 
et al. (2003) investigated the strength of the 6545 feature in 
V1016 Cyg, and suggested the opening half angle 6s ~ 60° 
adopting the distance estimate by Brocksopp et al. (2002). 

We adopt /i = 1.4 as suggested by Miirset & Nuss- 
baumer (1994), Voo = 10 km s~^. We choose a sightline from 
the white dwarf with 6o just a little bit smaller than 6s, for 
which p ~ d sin 6*3 ~ 73 AU. We may identify the H I col- 
umn density N%j = 1.2 x lO^'^ cm~^ that yields the observed 
centre shift with that for this sightline. We note that the ex- 
act choice of So affects the final result by a factor of two. 
However, in the case where the Stromgren boundary is well 
approximated by a hyperbola as noted by Schmid (1995), 
we believe that So is not much greater than p. Here, we just 
neglect this term and set Nhi — M /^fimpVooP]- Prom this 
consideration, we may obtain 

M = 3.6 X 10"'' Mq yr"\ (21) 

Because our analysis in the previous section is based 
upon a single column density medium, and each sight line 
toward the H I region as shown in Pig. 9 will be characterized 
by a higher H I column density than N%j. Therefore, the 
contributions to the centre shift from these sightlines will 
be smaller than that from the sightline for the Stromgren 
boundary. In order to account for the observed redward shift, 
we have to lower the column density, which will lead to a 
lower mass loss rate. However, this result should be taken 
with much care, noting that the estimate of Nhi is highly 
sensitive to the kinematics between the emission region and 
the neutral scattering region. 



5 DISCUSSION AND OBSERVATIONAL 
RAMIFICATIONS 

The inaccuracy of the estimation of Nhi can be attributed 
to a number of assumptions adopted in this work. The first 
assumption is that the scattering region is static with re- 
spect to the emission region. It is very plausible that the 
neutral region around the giant component takes the form 
of the slow stellar wind with a typical velocity of 10 km s~^. 
The amount of center shift 0.64 A in the Raman scattered 
He II A 4850 may also be obtained from the receding veloc- 
ity V = 7.9 km s^^ of the scattering region with respect to 
the emission region. This velocity scale is smaller than the 
slow stellar wind velocity, and hence our conclusion may be 
significantly affected. This is much more serious in the case 
of RR Tel, for which a smaller wavelength shift is concerned. 

However, it is also expected that the neutral region pos- 
sesses both receding and approaching motions that will com- 
pensate the kinematic effect to some degree. A more refined 



analysis will be possible if we can secure other Raman fea- 
tures including He II A 4340 formed blueward of H7, because 
the amounts of the centre shifts will differ if the effect is due 
to atomic physics. Currently, it is not feasible to determine 
the center location of this feature in our spectra. Spectra 
with longer exposures will be useful in establishing a consis- 
tent model. 

The assumption that the scattering region is character- 
ized by a single column density may not be valid. A more 
realistic picture would be a system consisting of spherical 
shells with density decreasing as r~^, where r is the dis- 
tance from the centre of the giant. However, the Stromgren 
boundary due to the strong UV radiation from the white 
dwarf region complicates this simple picture. Therefore, a 
more realistic model should be obtained by combining pho- 
toionization calculations that should also yield correct emis- 
sion line ratios (e.g. Miirset & Nussbaumer 1994). 

The measurement of the centre shift provides an upper 
limit of the mass loss rate with the assumption that the 
scattering region is static with respect to the He II emission 
region. On the other hand, Lee et al. (2003) provided a rough 
estimate of the lower bound M = 3 x lO"'^ M0 yr"'^ for 
the mass loss rate by equating Nhi of the sightline for the 
Stromgren boundary with the column density corresponding 
to the unit Raman optical depth. We tentatively conclude 
that the true mass loss rate may be found between the lower 
and upper bounds found by JL04 and the current work. 

When van Groningcn (1993) reported the discovery of 
the Raman scattered He II A 4850 feature in his spectrum of 
the symbiotic nova RR Tel, he also noted that it exhibits a 
similar redward centre shift by an amount Av — 12 km s""*^ 
in the parent velocity scale. He attributed this centre shift to 
a relatively receding motion between the emission region and 
the scattering region. He adopted the atomic data different 
from ours, where the line center of He II A 972 appears at 
972.09 A in his analysis. 

If we interpret his result using the atomic data adopted 
in this paper, the redward centre shift of the Raman scat- 
tered He II A 4850 in RR Tel amounts to AA = 0.33 A, a 
significantly smaller than the value AA ~ 1 A he quoted. 
If the scattering region in RR Tel is static with respect to 
the He II emission region, this redward shift indicates the 
H I column density Nhi — 4 x lO'^^ A, which is consider- 
ably higher than that in V1016 Cyg. This example shows 
that the analysis of Raman scattered features requires ex- 
tremely accurate atomic data and that Raman features can 
be a very sensitive probe of the density structure around the 
giant component in symbiotic stars. 
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